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ABSTRACT 

T/ST ACS/HRC images in UV (F250W), V (F555W), and I (F814W) resolve three isolated OB associations that 
lie up to 30 kpc from the stellar disk of the SO galaxy NGC 1533. Previous narrow-band Ha imaging and optical 
spectroscopy showed these objects as unresolved intergalactic H II regions having Ha luminosities consistent with 
single early-type O stars. These young stars lie in stripped H I gas with column densities ranging from 1.5 - 2.5 
xlO^" cm~^ and velocity dispersions near 30 km s~^. Using the HST broadband colors and magnitudes along with 
previously-determined Ha luminosities, we place limits on the masses and ages of each association, considering the 
importance of stochastic effects for faint (My > —8) stellar populations. The upper limits to their stellar masses 
range from 600 to 7000 M©, and ages range from 2-6 Myrs. This analysis includes an updated calculation of 
the conversion factor between the ionizing luminosity and the total number of main sequence O stars contained within 
an H II region. The photometric properties and sizes of the isolated associations and other objects in the HRC fields 
are consistent with those of Galactic stellar associations, open clusters and/or single O and B stars. We interpret 
the age-size sequence of associations and clustered field objects as an indication that these isolated associations are 
most likely rapidly dispersing. Furthermore, we consider the possibility that these isolated associations represent the 
first generation of stars in the H I ring surrounding NGC 1533. This work suggests star formation in the unique 
environment of a galaxy's outermost gaseous regions proceeds similarly to that within the Galactic disk and that star 
formation in tidal debris may be responsible for building up a younger halo component. 

Subject headings: galaxies: star clusters - H II regions - intergalactic medium - galaxies: halos - stars: 
formation 



1. INTRODUCTION 

H II regions in gaseous tidal arms, far beyond the ro- 
tating disks of galaxies, present an opportunity to study 
the star formation process in conditions vastly different 
from those in typical galactic spiral arms. They may 
even contribute to our understanding of satellite galaxy 
formation and evolution, halo stellar populations, and 
the enrichment of the intergalactic medium. Such inter- 
galactic H II regions have been found in galaxy clusters, 
compact groups, and in the halos of galaxies (e.g. Bo- 
quien et al. 2007; Walter et al. 2006; Ryan-Weber et 
al. 2004; Oosterloo et al. 2004; Mendes de Oliveira 
et al. 2004; Cortese et al. 2003; Sakai et al. 2002; 
Gerhard et al. H II 2002; Arnaboldi et al. 2002; Gal- 
lagher et al. 2001). In each case, they represent young 
stellar associations forming ostensibly where no stars ex- 
isted previously. Intergalactic H II regions are distinct 
from the H II regions in the outer arms of spiral galax- 
ies (e.g., Ferguson et al. 1998) in that they are further 
from the host galaxy, they do not show significant contin- 
uum emi ssion, and thus h ave much higher Ha equ i valent 
widths dRvan- Weber etr al. 200^; iGerhard et all 120021 : 
iMendes de Oliveira et al...2004) . These regions may rep- 
resent Ha-emitting counterparts to the extended UV- 
disk population of star clusters revealed by GALEX, 
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as they lie at similar distances from the host galaxies, 
and e xhibit similar levels of star-formation ('T hilker et al.l 
l2007t l. Much speculation exists as to whether intergalac- 
tic H II regions represent isolated stars, what role tidal 
interactions play in their formation, whether their stellar 
populations are typical of young star clusters in quiescent 
systems, and whether they contribute significantly to the 
enrichment of the IGM. 

Five intergalactic H II regions, three of which have 
been spectroscopic a,lly confirmed, lie in the vi cinity of the 
galaxy NGC 1533. iRvan- Weber et al.l (|2004D discovered 
these H II regions in the continuum-subtracted narrow- 
band image of NGC 1533 taken by the Survey for Ion- 
ization in Neutral Gas Galaxies (SINGG, Meurer et al. 
2006). The peculiar ring-like distribution of H I around 
NGC 1533 as seen in Figure 1 (and figure 2 of Ryan- 
Weber et al. 2004 with higher H I cut-off levels) may be 
the result of an interaction with the nearby dwarf galaxy 
seen in the NW corner, IC2038 (with the SINGG identi- 
fication J0409-56:S2). No obvious optical counterpart to 
the disturbed H I appears in either the DSS image or the 
SINGG image. The velocities of the three spectroscop- 
ically confirmed regions coincide well with the H I gas, 
which is bound to NGC 1533 (fRvan- Weber et al.ll200l . 
The projected radius from the center of NGC 1533 to 
the isolated H II regions ranges from 19 kpc (region 5) 
to 31 kpc (regions 1 and 2), up to four times the R-band 
25*'' magnitude isophotal radius of NGC 1533. Table 1 
summari zes the properties o f the t hree H II regions pre- 
sented in iRvan- Weber et al.l (|2004f ) and imaged here. It 
gives the regions' positions, projected radn from NGC 
1533, velocities determined from their Ha emission lines. 
Ha luminosities, and H I columns at each region's posi- 
tion. Note that the values in our Table 1 differ from those 
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Fig. 1. — The Ha SINGG image of NGC 1533 overlaid with Austraha Telescope Compact Array (ATCA) H I contours and marked with 
the locations of the individual spectroscopically confirmed H II regions studied here. The detection images are shown to the right of the 
SINGG image. The ATCA H I contours are 1.0 (in bold; 3a-), 1.5, 2.0, 2.5, and 3.0 xlO^^ cm"^ and have a resolution of ~ 1 arcminute. 
There is no detectable H I emission on the optical core of NGC 1533. 



presented in iRvan- Weber et al.1 (|2004[ ) because of the re- 
cent revision in the measurement of the distance to NGC 
1533 using images taken in paraUel to those discussed 
here (see Barber DeGraaff et al. 2007). In this paper, 
we present Hubble Space Telescope Advanced Camera for 
Surveys/ High Resolution Channel (ACS/HRG) images 
and photometry of the three spectroscopically-confirmed 
isolated H II regions associated with NGC 1533. These 
observations enable us to resolve the star-forming re- 
gions, identify their ionizing sources, and carry out an 
in-depth study of their stellar populations. 

In combination with the HST data, we make use 
of an ATCA H I syn thesis map of NGC 1533 
(jRvan- Weber et aDl2003bf ) to help constrain some of the 
details of the formation and evolution of the isolated H II 
regions. Intergalactic H II regions are often found form- 
ing in relatively low column density H I gas. There is 
no H I detected to a limit of 2 xlO^^ cm~^ at the lo- 
cation of the H II region near NGC 4388 (Oosterloo & 
van Gorkom, 2005), and the H I surface density in the 
vicinity of the three confirmed systems, regions num- 
bered 1, 2, and 5 in Fig ure 1 is 1.5 — 2.5 x 10^" cm~^ 
(jRvan- Weber et al]l2004l ). Furthermore, the H I peaks 
in the neutral gas ring surrounding NGC 1533 and in the 
stripped material around NGC 4388 do not exhibit any 
detectable level of star formation. These observations 
offer further support for the lack of correlation between 
star formation rate densities and H I surfac e densities on 
local scales found bv lKennicutt et al] (|2007f ). unlike what 
is seen on global galactic scales. The overall kinematics 
of the gas in the region of the intergalactic H II regions 
may provide clues as to the trigger of the star formation 



and the potential growth of the stellar associations into 
a tidal dwarf galaxy. 

Independent of the gaseous properties of cluster envi- 
ronments, star-cluster "infant mortality" suggests that 
the isolated, young star clusters prese nted here m ay be 
short-lived (iRafelski fc Zaritskvl ^OOSl : iFall et all [20051 : 
iChandar et all 120061 WiU supernova and mass loss in 
these clusters remove enough interstellar matter to leave 
the stars freely expanding, resulting in the clus ters' rapid 
disso lution? St udies of the Antennae galaxies ijFall et al.l 
I2005t fChandar ^eTaLl l2006l) and the Small Magellanic 
Cloud (Rafclski & Zaritsky"2005') purport that over half 
of clusters formed dissolve less than 10 Myr after their 
formation independent o f their initia l mas ses. How- 
ever, a recent study by iGieles et al.l (|2007f ) concludes 
that there is little or no "infant mortality" for ages be- 
yond 10-20 Myrs in the S mall Magellanic Cloud, and 
Ide Griis fc GoodwinI (|2007D find less than 30% of young 
(< 20 Myrs) SMC clusters undergo "infant mortality." 
Generally, studies of this kind require a statistical sam- 
ple of clusters. Instead, resolving a few star clusters in 
an isolated environment in a period during which they 
may be actively dissolving offers a different approach to 
examining their fates. 

The photometric properties of the underlying stars 
in these intergalactic H II regions can help to answer 
some of the questions surrounding their formation and 
evolution. Previous studies have used stellar popula- 
tion synthesis models and evolutionar y tracks to infer 
cluster ages and masses. For example, IWhitmore et al.l 
(1999) identify 3 distinct populations of star clusters in 
and around NGC 4038/4039 (Antennae Galaxy) using 
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TABLE 1 

Properties of the Confirmed Isolated H II Regions in NGC 1533 



H II RA Dec Separation Velocity Log L^f^ ^ N// / 

Region kpc km s^^ ergs s^^ lO^^cm"^ 

1 04 10 13.7 -56 11 37.5 31 846 ± 50 37.67 ± 0.26 2.4 

2 04 10 14.5 -56 11 35.9 31 831 ± 50 37.50 ± 0.26 2.4 
5 04 10 15.7 -56 06 16.2 19 901 ± 50 37.35 ± 0.26 1.5 



^ Summary of properties originally presented in Ryan- Weber et al. (2004): Positions are given 
in J2000 coordinates, with RA in hours, minutes, and seconds, and Declination in degrees, 
arcminutes, and arcseconds. Separation is the projected separation in kiloparsecs to the optical 
center of NGC 1533. The H I column densities are measured from ATCA maps, at the position 
of each H II region. NGC 1533 has a heliocentric velocity measured by HIPASS to be 785 km 
s~^. The numbers for projected separation and Ha luminosity differ from those presented in 
IRvan- Weber et al because we use the most recent measurement of the distance to NGC 

1533, 19.4 Mpc (Blakeslee, private communication).^ As measured in the SINGG Ha images 



their integrated photometric properties. iGerhard et al.l 
()2002t l determine the mass and age of an isolated H II 
region in the Virgo Cluster from its measured Ha and 
V band luminosities. Thus far, studies of these sorts 
of objects have focused on what the star clusters' inte- 
grated light reveals about their formation histories, and 
have been unable to resolve the clusters into individual 
stars or stellar components. With the discovery of more 
and more young, faint clusters in tidal tails of merging 
galaxies, intracluster space and the intergalactic medium, 
there arises a need for information about the nature of 
the forming stars. Yet, there is an unavoidable, inher- 
ent uncertainty in the properties of lower mass clusters 
derived from models that a ssume a fully populated IMF 
(jCervino &: Luridianall20d^ . Resolving the star clusters 
presents an opportunity to examine their populations in 
greater detail. 

We organize this paper as follows: Section 2 contains a 
description of our observations and measurements; Sec- 
tion 3 presents the results and model comparisons; Sec- 
tion 4.1 explores comparisons between these intergalactic 
H II regions and other extragalactic and Galactic young 
stellar populations; in Sections 4.2 and 4.3, we discuss 
the gas properties in the vicinity of the intergalactic H II 
regions and examine the likely fates of the stellar associ- 
ations; and Section 5 briefly summarizes our results. The 
Appendix contains the new calculation of the conversion 
factor between the ionizing luminosity and the total num- 
ber of main sequence O stars contained within an H II 
region following the methodology of Vacca (1994). 

2. OBSERVATIONS AND MEASUREMENTS 

2.1. HST Observations 

HST ACS/ HRC observations were carried out in Oc- 
tober 2004. The targets were observed in three filters: 
F250W (5808 s); F555W (2860 s); and F814W (2892 s), 
covering a wide wavelength range chosen to sample the 
shape of the ultraviolet to optical SED. Throughout this 
paper, we will refer to the F250W filter as UV, F555W as 
V, and F814W as I in the interest of brevity and simplic- 
ity. Each orbit was split into two exposures separated by 
a small dither, resulting in 5a limiting AB magnitudes of 
28.0 in UV, 28.9 in V, and 29.0 in I. The images were pro- 
cessed with the CALACS calibration pipeline in order to 
remove the instrumental signa ture, and then combined 
using the Apsis pipeline dBlakeslee et al.|[2003() resulting 
in cleaned and combined single images in each filter at 



each pointing at the optimal pixel scale of ACS/HRC, 
0'.'025 pixel" ^. We create a detection image by sum- 
ming the images in the three filters. Due to the paucity 
of known point sources in our HRC images, we mea- 
sured the PSF from HRC images taken in the same 
cycle and in the same filters of the center of 47 Tuc, 
using at least 50 stars per image. The point sources 
in each image have an average Gaussian FWHM of 2.6 
pixels, independent of the 3 filters investigated. At the 
adopted distance to NGC 153 3, assumed to be 19.4 Mpc 
(|Barber DeGraaff et ani2007f) . this FWHM corresponds 
to a physical size of 6.1 pc. 

Figure 1 shows the ACS/HRC detection images of the 
OB associations powering the isolated H II regions along 
with their locations i n the R-band cont inuum SINGG 
image of NGC 15 33 (iMeurer et al.ll200l); w e show the 
H I distribution (|Rvan- Weber et al.l l2003bf ) with con- 
tours. The half-light radii of associations 1, 2, and 5 
are 7.9, 9.0, and 6.2 pixels respectively, corresponding to 
physical sizes of 18.5, 21.1, and 14.5 pc. At a distance of 
19.4 Mpc, 1 arcsecond corresponds to 94 parsecs. Three 
to nine clumps, most close in size to the image PSF, com- 
prise each association. While multiple sources bright in 
UV and V surround associations 1 and 2, no such objects 
appear within the same proximity to association 5. 

2.2. Source Detection and Definitions 

We catalogued the sour ces in our im ages by applying 
the SExtractor code (Berti n fc Arnout s 1996), as imple- 
mented in the Apsis pipeline, to the relevant detection 
images. We eliminated spurious sources (e.g. single pixel 
detections, image edge artifacts) and objects detected 
in only one band by inspecting the SExtractor detec- 
tions by-eye. Because the bulk of the following analysis 
requires color measurements, single-band detections (4 
sources) are not of particular usefulness to us, though 
they may indeed be real. Of the initial 30 sources iden- 
tified by SExtractor in the image containing associations 
1 and 2, we include 19 in this analysis. Of the initial 
36 sources identified by SExtractor in the image contain- 
ing association 5, 11 appeared to be genuine detections, 
including one single-band detection (F814W) of an ob- 
vious background galaxy. Figure 1 shows the HRC de- 
tection images containing the labelled isolated associa- 
tions. SExtractor split association 1 into three separate 
objects, and detected associations 2 and 5 as single ob- 
jects. Further examination revealed each association to 
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Fig. 2. — Cut-out images in each HRC filter of the three isolated associations and their components (which wc refer to as "clumps"), 
circled and labeled. The image scale is shown in the top left corner. F250W is a UV bandpass, while F555W and F814W closely correspond 
to V and I bands, respectively. We present a detailed description of source detection and definition in Section 2.2. 



be composed of multiple chimps insufficiently deblended 
by SExtractor. Within the associations, we separated 
the components using IRAF contour plots with levels of 
2(7, 4(7, 6(7 and 8(t, defining an individual component as 
having at least two contours, corresponding to at least 
4(j above; the background. These criteria resulted in six 
components in association 1, nine components in associ- 
ation 2, and three components in association 5 (see Fig- 
ure 2). Finally, using the IRAF task IMEXAM, we fitted 
Gaussian distributions to all field objects and clumps in 
order to determine their FWHMs. 

There are three main types of objects to which we will 
refer throughout this paper: the stellar associations pow- 
ering the total H II regions (=associations), the compo- 
nents of the associations (=clumps) , and the field sources 
in the vicinity of the associations (= field objects). Fur- 
thermore, we assume that all of the extended and red 
field objects in the vicinity of the associations are back- 
ground galaxies (see Section 2.4 and Table 4), and thus 



exclude them from most of this analysis. We address 
background galaxy contamination of the rest of the field 
objects in Section 2.4. The images of the objects included 
in this study are presented in Figures 1-4. Figure 1 gives 
an idea of the environment containing the intcrgalactic 
associations. All three lie in a low column density H I 
ring, with projected radii 19 - 31 kpc from the nucleus of 
NGC 1533. Figure 2 displays close-ups of the clumps in 
the three bands, UV (left column), V (center column), 
and I (right column), circled and labeled. The letters in 
the images correspond to the names in Table 2. Figure 3 
shows the two HRC fields (in the three combined filters), 
with labels next to each object. Finally, Figure 4 shows 
three-color image close-ups of the associations, clumps, 
and field objects in the two sets of HRC images. In these 
color close-ups, blue, green, and red represent the UV, 
V, and / bands, respectively. 

2.3. Aperture Photometry 
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Fig. 3. — HRC detection images show the full fields containing the isolated stellar associations 1, 2, and 5. The 30 detected field objects 
are labeled according to their names in Tables 2, 3 and 4. 
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Fig. 4. — Three-color images of associations and field objects where blue, green and red represent UV (F250W), V (F555W), and I 
(F814W), respectively. The objects are labeled with names corresponding to those given in Tables 2, 3, and 4. 
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TABLE 2 

Photometry Results and Sizes for Isolated Associations and Their 

Components 



Assn 




r A (i^ zaUW } 
10" cgs'' 


■\ /■ Ta 
V-i 

Vega 


UV-V 
AB 


my 




Mv 


r WrliVl 
pc 


1 


Total 


8.13 ±0.78 


-o.i2t;lE! 


-0.211:1^ 


23.381 


.04 
.04 


-8.06 


61 




a 


0.24 ± 0.04 




-0.861-30 


27.851 


.17 
.20 


-3.59 


3.39 




b 


0.52 ±0.07 


-0.22t;2^ 


-0.481-11 


26.641 


.09 
.10 


-4.80 


4.91 




c 


0.18 ±0.04 


-0.33+-3g 


-o-isl-JI 


27.5ll 


.13 
.15 


-3.93 


5.48 




d 


2.17 ±0.14 




0.331;°^ 


24.271 


.02 
.02 


-7.17 


9.40 




e 


0.21 ±0.04 


-0.371 


0.07+-2? 


27.061 


.09 
.10 


-4.38 


4.00 




f 


1.95 ±0.09 


-0.451:1° 


-0.66i;oe 


25.381 


.04 
.04 


-6.06 


4.98 




DLC 


2.63 ±0.78 


-0.25til 


0.33l:i 


24.041 


.09 
.10 


-7.40 




2 


Total 


11.2 ±0.77 


0.08+-fo 


-0.67l;;j« 


23.511 


.04 
.05 


-7.99 


67 



a 


0.26 ± 0.04 


-0.02+-27 


-0.541 


.21 
.24 


27.471 


.12 

.14 


-3.97 


4.44 


b 

c 


1.18 ±0.07 
1.40 ±0.07 


-0.31+-;^ 


-0.471 
-0.671 


.07 
.08 
.06 
.07 


25.731 
25.751 


.04 
.04 
.04 
.04 


-5.71 
-5.69 


9.66 
8.67 


d 


0.62 ± 0.05 




-0.831 


.11 
.12 


26.80l 


.07 
.08 


-4.64 


5.31 


e 


1.25 ±0.07 


-0.25+-?9 


-0.621 


.07 
.08 


25.821 


.04 
.05 


-5.62 


9.60 


f 


0.93 ± 0.05 


-0.45l-:i^ 


-0.741 


.08 
.08 


26.271 


.05 
.05 


-5.17 


< 4.70 


g 


0.85 ±0.07 


-0.171-lJ 


-0.641 


.10 
.11 


26.261 


.06 
.07 


-5.18 


4.83 


h 


0.17 ±0.04 


-0.091-JI 


-O.lll 


.27 
.34 


27.471 


.12 
.13 


-3.97 


< 4.70 


i 


0.58 ±0.07 


0.19l|0 


-i.osl 


.17 
.18 


27.061 


.12 
.13 


-4.38 


5.03 


DLC 


3.96 ± 0.77 


0.171:1^ 


-0.261 


.20 
.24 


24.191 


.10 
.11 


-7.25 




5 Total 


3.51 ±0.66 


-0.221;!^ 


0.021 


.20 
.23 


24.061 


.06 
.06 


-7.38 


40 



a 


0.15 ±0.05 


0.46l:^« 


-0.55l-i^ 


28.051 


.20 
.25 


-3.39 


4.46 


b 


0.61 ±0.07 


-0.45+-!^ 


0.021-^^ 


25.961 


.05 
.06 


-5.48 


< 4.70 


c 


1.57 ±0.10 


-0.54l-i? 


-o.33i-.E;| 


25.281 


.04 
.04 


-6.16 


9.54 


DLC 


1.33 ±0.66 


-0.541-J5 


0.74l:^,^ 


24.371 


.10 
.11 


-7.07 





^ V - I colors have been transformed to Vega Magnitudes in the standard Johnson-Cousins 
filter system using the prescription of Sirianni et al. (2005). V band magnitudes are 
technically in units of Vega magnitudes, although magnitude systems converge in the V 
band, making the given values roughly correct for the AB magnitude system as well. The 
UV - V colors are left in AB magnitudes, as there is no calibration for the conversion 
of the F250W filter to Vega magnitudes. Absolute magnitudes were calculated using a 
distance modulus of 31.44 (Blakeslec, private communication). For each association, V 
band magnitudes were corrected for nebular emission ([OIII]AA4959,5007).'^ ergs s""^ cm~^ 



We performed circular aperture photometry on each 
source using the IDL routine APER. Aperture sizes for 
each association and field object "were determined by ex- 
amining where the curve of growth (in flux vs. aperture 
radius) began to flatten, whereas aperture sizes for the 
clumps were determined by eye. For each association and 
field object, sky annuli were positioned approximately 1 
arcsecond away from each of the objects' centers. Be- 
cause the clumps are so closely spaced in the image, we 
subtracted the sky measurement outside the correspond- 
ing association as a best estimate. As evidenced by Fig- 
ure 2, the circular apertures containing the clumps do 
not account for all the light in the association. We chose 
the aperture sizes to contain the maximum amount of 
emission while not overlapping with other apertures. By 
masking out the apertures containing the clumps, and 
running APER on what remained, we determined fiuxes 
for the diffuse light components (DLC) of each associa- 
tion. This component simply represents the light losses 
resultant from the small apertures, and thus may not be 
genuine diffuse emission. Thus, the DLC is essentially 
the difference between the total association fiux, and the 
sum of the fiuxes of the individual clumps. 



A correction for f oreground extinction, based on the 
ISchlegel et alj (|1998D extinction maps, (foreground E(B- 
V) = 0. 016), w as applied according to the parameteri- 
zation of iCarde lli ct al. (1989), including the update for 
the near-UV given bv lO'Donnelll (jl994l ). Given the loca- 
tion of the isolated associations far outside the optically 
luminous area of NGC 1533, we excluded any correction 
for the internal extinction. Additionally, fiuxes for the 
associations were corrected for the strong emission lines 
[O III]AA4959, 5007 present in the F555W bandpass. We 
used the STSDAS package SYNPHOT to simulate pho- 
tometry using an actual spectrum of association 1 (to be 
presented in a future paper) both with and without the 
nebular emission lines. The procedure for this correc- 
tion is as follows: (1) We measured the Ha fiux for all 
three H II regions in the NGC1533 SINGG narrow-band, 
continuum-subtracted image; (2) We measured the Ha 
emission line fiux in the spectrum of association 1; (3) 
We divide quantity (1) by quantity (2) for each associa- 
tion to obtain a flux correction factor for each individual 
association; (4) Using SYNPHOT, we measure the total 
flux of association I's spectrum contained in the F555W 
bandpass; (5) We subtract the emission lines in this area 
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of spectrum, and remeasure the total flux in the same 
area with SYNPHOT; (6) We subtract (5) from (4) to 
get an uncahbrated flux of just the emission hnes in the 
F555W bandpass; (7) Finally, we calibrate and scale the 
total nebular emission-line flux from (6) by the individual 
correction factor determined in (3). 

These emission-line corrections were 0.26, 0.17, and 
0.27 AB magnitudes (17% - 28% of the total flux) for 
associations 1, 2 and 5, respectively. We note that these 
emission line corrections assume that the flux ratios of 
Ha to [O III] emission lines are the same for each associ- 
ation. The error due to potential reasonable variations in 
these line ratios is negligible when compared to the total 
corrections. Furthermore, this correction has the effect 
that the associations are redder in V-I after correction 
than one would expect based on the V-I colors of the 
individual clumps, which were not corrected for nebular 
emission due to the limited (~1.5 arcsecond) resolution 
of the SINGG Ha images. 

Finally, for comparison with population synthesis mod- 
els and previous work, we transformed the ST magni- 
tudes to the Vega magnitude system, and the filters 
to the standard Joh nson-Cousins system following the 
method outlined by ISirianni et al.l (2005) for both the 
F555W and F814W data. Because magnitude systems 
converge in the V band, and because the F555W filter is 
quite similar to the Johnson-Cousins V band, its trans- 
formation produces virtually no change. Furthermore, 
the F250W filter has no Johnson-Cousins counterpart, 
and it is impractical to convert AB magnitudes to Vega 
magnitudes in this part of the spectrum. Throughout 
this paper, we will quote V-I colors in terms of Vega mag- 
nitudes in the standard Johnson-Cousins filters. This 
magnitude system allows for the greatest ease in com- 
paring our results with both previous work and existing 
population synthesis models. However, UV-V and UV-I 
colors and UV magnitudes will remain in the HRC fil- 
ter bandpasses and in AB magnitudes. The faintness of 
the clumps and field objects combined with considerable 
pixel-pixel noise in the read noise limit produces errors 
in the photometry as high as 0.25 magnitudes. The re- 
sulting Vega magnitudes are presented in Tables 2, 3, 
and 4 as well as the V-I and UV-V colors, along with the 
measured FWHM of each source. 

2.4. Source Properties and Contamination by 
Background Galaxies 

Our sample of field objects and clumps totals 48 indi- 
vidual sources altogether, 34 of which are located in the 
3 kpc X 3 kpc field around associations 1 and 2. We 
estimate that the associations span 40 - 75 parsecs in 
their full visible extent, similar to large Milky Way stel- 
lar associations. Figure 3 gives the location of all objects 
identified in the two HRC images, with the field contain- 
ing associations 1 and 2 on the top frame, and that of 
association 5 on the bottom frame. Of the 48 individual 
sources, only 17 are resolved. 9 of the 11 field objects in 
the association 5 field are well-resolved extended sources. 
The V-I colors of the clumps range from -0.7 to 0.2 in 
V-I, and those of the field sources range from -0.9 to 1.8. 
If all the sources were at the same distance as NGC 1533, 
the absolute y-band magnitudes range from -3.4 to -7.3. 

Almost certainly, some of the field objects in the HRC 
images represent distant background galaxies. From 



Hubble Deep Field galaxy number densities, we ex- 
pect 1-2 X 10^ galaxies per square degree in WFPC2, 
corr esponding to a limiti ng I-band magnitude of 25- 
26 (jCasertano et al.ll2000f ). Field galaxy number den- 
sities are lower in the far-UV, roughly 3 x 10"* galaxies 
per square degree t o a limiting AB magnitude of 26.5 
^Teplitz et al.ll200"6l ). Considering these number densi- 
ties, we would then expect 5-10 red background galaxies 
to appear in our 29"x26" field of view, and 0-1 galax- 
ies with emission in the UV. We detect and analyze 11 
objects in the vicinity of association 5, none of which 
appear in the UV images, and the majority of which 
appear quite red in V-I. Of the 11 objects in the vicin- 
ity of association 5, six are very obviously extended or 
elliptical red sources, with V-I colors (Vega magnitudes) 
between 0.7 and 1.8 (as previously mentioned, one object 
is detected only in F814W). Three of the remaining five 
compact sources have V-I colors in a similarly red range, 
while two exhibit fairly blue V-I colors, 0.04 and 0.19. 
The lack of UV emission from these sources combined 
with the number being consistent with our prediction for 
background galaxy detections leads us to classify all 11 
as background emitters. 

We detect and analyze 19 objects in the vicinity of as- 
sociations 1 and 2. In contrast to the association 5 field 
objects, many of these field objects appear bright in UV 
and V, suggesting that the star formation in the SE part 
of NGC 1533's H I ring extends beyond the initially de- 
tected isolated H II regions. Furthermore, most of the 
surrounding blue field objects lie clustered within 500 
parsecs of the nearest association (1 or 2), while the likely 
background galaxies are spread more evenly throughout 
the image. While we cannot say for certain which of the 
field objects around associations 1 and 2 are background 
galaxies, we can predict that the 10 sources not detected 
in F250W, 5 of which have V-I colors greater than 1.0, 
are more likely to be distant background galaxies. These 
sources also appear to be more extended. Moreover, this 
number of sources is consistent with the prediction based 
on the background galaxy density of deep WFPC2 im- 
ages. 

We do note, however, that we are unable to distinguish 
extended likely background galaxies from potentially old 
(~ 1 Gyr), extended open clusters, which can have V-I 
colors up to 1.4 and abs olute ma gnitudes ranging from 
-1 to -5.1 in the V band ()Lata e t al. 2002). Nor can we 
differentiate unresolved likely background galaxies from 
globular clusters, which have V-I colors ranging from 0.85 
to 2 65, and abs olute F-band magnitudes from -1.5 to - 
9.5 ()Harri£lll996l ). Still, such old, extended open clusters 
are rare in the Milky Way (e.g. Lata et al. 2002; Lada 
& Lada 2003), and globular clusters tend to lie closer to 
the galaxy center (within^ 20 kpc) than the star for- 
mation seen here ( Harris! [1996( 1. although a few are as 
far as 100 kpc from Galactic center. Still another possi- 
bility remains that some of these objects could represent 
isolated red supergiants. 

Tables 3 and 4 present the photometric properties and 
optical sizes for all detected field objects. In these tables, 
we have identified those field objects which are the best 
candidates for distant background galaxies (extended, 
red, and/or no UV-emission) as F, and those blue field 
objects which are most likely to be associated with the 
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ongoing star-formation around NGC 1533 as F. Figure 
3 shows three-color (RGB) images of the field objects 
and background galaxies. The top panel shows the field 
objects surrounding associations 1 and 2, and the bot- 
tom panel displays the field objects of association 5. The 
labels on each image correspond to the names listed in 
Tables 3 and 4 ("Assn" for the association; "F" for field 
object, and "F" for likely background galaxy). We do 
include all of these objects in the subsequent analysis for 
completeness' sake (except where their red colors push 
them outside of adopted plotting ranges), and we urge 
the reader to consider the 10 "dotted" objects in Table 
3 and the 11 objects in Table 4 (also identified in the 
figures) as potential background galaxies. 

In Figures 5 and 6, we compare the individual sources' 
absolute magnitudes and V-I and UV-V colors with 
those of single star MK spectral types and a 50 Mq 
instantaneous-burst Starburst99 population (see 3.1). In 
order to plot MK spectral types in My versus UV-V, we 
converted their U-band absolute Vega magnitudes to the 
F250W filter in AB magnitudes using the SYNPHOT. 
To arrive at Starburst99 UV-V colors we calculated the 
F250W AB magnitude using the output synthetic spec- 
trum along with SYNPHOT. In many cases, the colors 
and absolute magnitudes of the objects are well-matched 
by single O or B stars. This similarity does not prove that 
field objects and clumps are single stars, but that O and 
B stars dominate the radiation output in many cases. In 
other cases, the expected magnitude-color evolution of a 
Starburst99 population (run for an instantaneous burst 
Salpeter IMF witha = 2.35, M„p = 100, Mion, = 0.5, 
Padova group isochrones, and Z= OAZq) with an initial 
mass of 50 Mq matches the data. We note that a 50 
Mq population with a 100 Mq upper mass limit is un- 
physical, and results in fractional stars. We have simply 
used an initial starburst mass of 50 Mq to appropriately 
scale the absolute magnitude and illustrate a degener- 
acy. The colors are unaffected by the mass of the initial 
starburst. We discuss the implications of fractional stars 
and stochastic IMFs in section 13.31 For completeness, 
we have included the bluest (in V-I) background galaxy 
candidates in Figure 5, using NGC 1533's distance mod- 
ulus to calculate their absolute magnitudes. The reddest 
candidates do not lie within the ascribed color range. 

3. STELLAR POPULATIONS OF ASSOCIATIONS 

Ideally, stellar population synthesis models, along with 
integrated colors and luminosities, can provide a solu- 
tion for the masses and ages of the stars populating 
the clusters. These models aim to account for most of 
the known stellar evolutionary processes, and integrate 
over various sets of assumed initial conditions. In the 
following analysis, we describe two approaches we took 
in an effort to constrain the stellar populations of the 
isolated associations. We compare the associations' col- 
ors and luminosities t o output from Starburst99 models 
(|Leitherer et all 119990 . and we perform a simple, inde- 
pendent calculation using the associations' total UV and 
Ha luminosities. Both approaches employ single-age, 
single-metallicity models with a Salpeter IMF. While the 
Starburst99 model accounts for details of stellar evolu- 
tion, it cannot account for anything but a scaled-down 
IMF over the full range of stellar initial masses. Our in- 
dependent calculation is more quantized with regard to 
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Fig. 5. — Color-magnitude diagram for the blue field objects 
and clumps, with V-I in Vega magnitudes. The dashed and solid 
lines shown are calibrations of MK spectral types from Schmidt- 
Kaler (1982). The dotted line represents the V-I color evolution of 
a Starburst99 50 Mq population with Z= 0.4Zq. The population 
synthesis model was run for an instantaneous burst Salpeter IMF 
{a = 2.35, Mup = 100, Mi^^ = 0.5), and employs Padova group 
isochrones. 
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Fig. 6. — Color-magnitude diagram for the field objects and 
clumps, with UV-V in AB magnitudes. The dashed and solid lines 
shown are calibrations of MK spectral types from Schmidt-Kaler 
(1982). The dotted line represents the UV-V color evolution of the 
same Starburst99 population shown in Figure 5. 

the stellar populations, but does not account for some 
of the details of stellar evolution. In section 13.31 we dis- 
cuss the limited applicability of both models to our data 
given the low luminosities of the isolated associations, 
and masses below what is required for a fully-populated 
IMF. 

3.1. Comparison with Starburst99 Model 

To accurately designate a metallicity for the Star- 
burst99 model, we used an optical spectrum of associ- 
ation 1, estimating its metallicity to be roughly ~ OAZq 
(Werk et al., in preparation). Figure 7 plots various di- 
agnostics from ages 1-16 Myr for a Starburst99 syn- 
thetic population for Z= 0.4Zq, an instantaneous burst 
of IO^Mq, and a Salpeter IMF with a lower limit of 0.5 
Mq and an upper limit of 100 Mq. It also shows the 
corresponding values calculated for associations 1, 2 and 
5. The error bars indicate photometric errors. 

In Figure 7a, we show the temporal evolution of 
the ratio Log {Fho/ Fuv) for the Starburst99 model 
and the three associations. Ha and UV l uminosities 
serve as reasonable star formation indicators (jKennicuttj 
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TABLE 3 

Photometry Results and Sizes for Association 1 & 2 Field Objects 



ID 


R.A. 


Dec. 


Fx (F250W) 


V-P 


UV-V 


my 


Mv FWHM 




h m s 


d m s 


10-^^ cgs^ 


Vega 


AB 




pc 



Fl 


04 


10 


14. 


,54 


-56 


11 


44. 


.84 






- 






F2 


04 


10 


15. 


,51 


-56 


11 


43, 


,35 






- 






F3 


04 


10 


12. 


,73 


-56 


11 


45, 


.58 


1. 


,36 


± 


0. 


14 


F4 


04 


10 


14. 


,64 


-56 


11 


39, 


.58 






- 






F5 


04 


10 


15. 


10 


-56 


11 


37, 


.89 


1. 


,02 


± 


0. 


13 


F6 


04 


10 


14. 


,82 


-56 


11 


36, 


.00 


0. 


,90 


± 


0. 


14 


F7 


04 


10 


14. 


,24 


-56 


11 


36, 


.35 






- 






F8 


04 


10 


14. 


,80 


-56 


11 


35, 


.18 


0. 


,56 


± 


0. 


11 


F9 


04 


10 


14. 


,89 


-56 


11 


34. 


.69 


0. 


,75 


± 


0. 


13 


FIO 


04 


10 


13. 


,55 


-56 


11 


36, 


,42 


1. 


,59 


± 


0. 


14 


Fll 


04 


10 


14. 


,74 


-56 


11 


34. 


.77 


0. 


,35 


± 


0. 


14 


F12 


04 


10 


14. 


,41 


-56 


11 


33, 


.89 












F13 


04 


10 


13. 


,74 


-56 


11 


34, 


.02 












F14 


04 


10 


13. 


,70 


-56 


11 


33, 


.91 


0. 


,68 


± 


0. 


17 


F15 


04 


10 


14. 


,05 


-56 


11 


32. 


.99 












F16 


04 


10 


12. 


,86 


-56 


11 


33, 


.76 












F17 


04 


10 


13. 


,49 


-56 


11 


31, 


.89 












F18 


04 


10 


14. 


19 


-56 


11 


30, 


.99 


2. 


,41 


± 


0. 


16 


F19 


04 


10 


13. 


,85 


-56 


11 


27, 


,36 













0.05: 
1.38: 

-0.23: 

i.oi: 
-o.4o: 
-o.6i: 

2.03: 

0.34: 
-0.27: 

-o.io: 

0.79: 

1.49: 

-0.89: 

1.36: 

0.67: 
-0.32: 
-0.27: 

-0.01 



.23 
.26 
.32 
.43 

-.04 

-.04 
.17 
.19 

-.34 
.45 
.36 

-.51 
.10 

.11 



.25 
.18 
.21 
.07 
.08 
.24 
.30 
.27 
.34 
.52 
.93 
.27 
.34 
.29 
.36 
.20 
.24 
.08 
.09 
+ .27 
33 



1.17: 

-0.95: 
-0.74: 

-1.07: 

-0.67: 
-0.86: 
2.04" 



-0.44' 



-0.19' 



.19 
.23 

-.20 
.23 
-.22 
-.25 



-.25 
-.30 

.23 
-.26 

.13 
-.14 

.36 
-.55 



.31 

.39 



-.08 
.09 



25.52 
26.63 
23.92 
26.55 
26.35 
26.27 
25.24 
27.12 
26.41 
25.78 
24.56 
26.95 
26.94 
26.28 
27.28 
26.61 
25.77 
24.66 
26.92 



+ .14 
.16 
+ .28 
.39 
.03 
.03 
.14 
.16 
.15 
.17 
.15 
.17 
.10 
.11 
.17 
.20 
.15 
-.17 
.09 
-.10 
.04 
-.04 
-.20 
.25 
-.25 
.32 
.19 
.23 
.25 
.32 
.24 
.30 
-.09 
-.10 
-.04 
-.04 
.16 
.19 



■5.92 


11.38 


■4.81 


7.81 


■7.52 


5.91 


■4.89 


7.88 


■5.09 


5.60 


■5.17 


5.95 


■6.20 


6.02 


■4.32 


3.97 


■5.03 


5.62 


■5.66 


7.15 


■6.88 


5.48 


■4.49 


5.34 


■4.50 


5.64 


■5.16 


4.66 


■4.16 


5.60 


■4.83 


8.02 


■5.67 


5.19 


■6.78 


5.62 


■4.52 


5.78 



^ V-I colors have been transformed Vega Magnitudes in the standard Johnson-Cousins filter system using the 
prescription of Sirianni et al. (2005). Absolute magnitudes and physical sizes were calculated using a distance 
modulus of 31.44 (D — 19.4 Mpc, Blakcslee, private communication). ^ ergs s"""" cm~^ A""'" 



TABLE 4 

Photometry Results and Sizes for Association 5 Field 
Objects 



ID 


R.A. 


Dec. 


V-I 


m„ 


FWHM^' 




h m s 


d m s 


Vega 




pc 



F20 
F21 
F22b 


04 
04 
04 


10 
10 
10 


16. 
16. 
16. 


,52 
,86 
,00 


-56 
-56 
-56 


06 
06 
06 


24. 
04. 
11. 


.07 
.30 
.62 


1, 
1, 


.22+ 
.711 


.28 
.36 
.22 
.28 


25. 
25. 


,27l 
.641 


.25 
.33 
.21 
.27 


14. 
15. 
17. 


.68 
.47 
.48 


F23 


04 


10 


15. 


,06 


-56 


06 


18. 


.80 


0, 


.04+ 


.16 
.18 


26. 


■ OOi 


.09 
.10 


6. 


.95 


F24 


04 


10 


14. 


,77 


-56 


06 


17. 


,57 


0, 


.19t 


.16 
.18 


25. 


,99l 


.11 
.12 


9, 


.99 


F25 


04 


10 


14. 


,48 


-56 


06 


20. 


.22 


1, 


.56+ 


.21 
.25 


24. 


.30l 


.19 

.23 


24. 


.73 


F26 


04 


10 


14. 


,66 


-56 


06 


10. 


,95 


0, 


.85+ 


.07 
.09 


24. 


,80+ 


.06 

.07 


9, 


.91 


F27 


04 


10 


14. 


,04 


-56 


06 


23. 


.07 


2, 


,25+ 


.28 
.38 


27. 


,14+ 


.27 

.37 


5. 


.89 


F28 


04 


10 


14. 


,20 


-56 


06 


18, 


,31 


0, 


.7lt 


.11 
.11 


24, 


.811 


.08 

.08 


17. 


.90 


F29 


04 


10 


14. 


,02 


-56 


06 


16. 


.17 


1, 


.34t 


.13 
.13 


26. 


.091 


.11 
.12 


6. 


.36 


F30 


04 


10 


14. 


,96 


-56 


06 


00. 


,15 


1, 


.791 


.32 
.43 


24. 


,83l 


.30 

.42 


24. 


.97 



^ Physical sizes arc given in parsecs for straightforward comparison with 
clumps and other field objects. ^ We do not detect F22 in any filter other 
than F814W, although we consider it a real source, with an / band mag- 
nitude of 23.14 (Vega magnitudes). 



I1998D . While Ha traces the ionizing radiation from stars 
more massive than ISM©, the UV luminosity originates 
from stars more massive than 3Mq. The ratio Log 
[Fhol/ Fuv) decreases steeply with time, as massive ion- 
izing stars die out, and late- type O stars and B stars 
produce the majority of the UV flux. The best-match 
ages to the evolution of Log (FHa/Fuv) are 3.2 Myrs, 
4.5 Myrs, and 3.0 Myrs for associations 1, 2, and 5 re- 
spectively. In this plot, the horizontal error bars repre- 
sent the minimum possible error given the photometric 
errors. 



The evolution of UV-V, UV-I, and V-I colors of the 
same Starburst99 population is plotted in three color- 
color diagrams for the same range in ages in Figures 
7b, 7c, and 7d. The Starburst99 UV magnitudes in 
the HRC F250W bandpass were generated by the SYN- 
PHOT package for the output Starburst99 spectra. Seen 
in 7b, 7c, and 7d, associations 1, 2, and 5 are consistently 
bluer than the model colors, except for association 2 in 
V-I. Thus, we cannot constrain the ages using simply the 
colors of the associations. As we will discuss, the inability 
of the Starburst99 model to match the association colors 
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Fig. 7. — Evolution with time of the Ha to UV ratio and relevant optical color-color diagrams for the Starburst99 stellar population 
synthesis model at Z= 0.4Zq. The model was run for an instantaneous burst Salpeter IMF {a = 2.35, Af„p = 100), and employs Padova 
group isochrones. a) The Ha to UV ratio of each association with a symbol at the intersection point of the Starburst99 model. The inferred 
ages are 3.2 Myr, 4.5 Myr, and 3.0 Myr for associations 1, 2, and 5 respectively. The errors in the best-match ages correspond to the 
uncertainty due to the photometric errors alone, and do not include stochastic effects or model errors, b) UV-I versus UV-V; the dotted 
line shows Starburst99 values for ages between 2-10 Myrs c) UV-V versus V-I d) UV-I versus V-I 



most likely arises from small number sampling in which a 
few early-type stars dominate the output radiation. An 
incompletely-sampled IM F results in a signi ficant scatter 
of colors at a given age (jChiosi et al.lll988^ ■ and would 
explain the discrepancies seen in figures 7b, c, and d. 

From the Starburst99 model, we can also derive mass 
estimates based on the absolute V band magnitudes and 
ages of the associations. At a particular age, a 10^ Mq 
Starburst99 stellar population has a predicted absolute 
V band magnitude, Mv,S99- We then use a simple scaling 
relation between mass and absolute magnitude to derive 
the masses of the associations: M = io6+0-4(Mv,s99-Mv)^ 
where M is the mass of the stellar population and My is 
the absolute magnitude of the association given in Table 
2. The total stellar mass of each association is 1490, 1850, 
and 840 M© for associations 1, 2 and 5 respectively. We 
repeat this calculation for the I-band magnitudes. The 
total stellar mass from the I-band magnitude of each as- 
sociations is 1240, 2030, and 610 M© for associations 1, 
2 and 5, very similar to the results of the V-band cal- 
culation. We note that this particular way of estimating 
the mass assumes a scaled-down IMF containing the full- 
range of stellar masses from 0.5 M© to 100 M©. Indeed, 
if we were to try to quantize the numbers of stars by their 



spectral type we would find only fractional stars for the 
highest masses. Table 5 includes columns for the ages 
and masses derived using the Starburst99 model. 

Although these plots and derived parameters are in- 
structive, they are limited by the low luminosities of the 
clusters. A very basic limit for the application of evolu- 
tionary synthesis data to observational data, termed the 
"lowest luminosity limit" (LLL) by Cervifio & Luridian^ 
(|2004l ). is that the total luminosity of the cluster mod- 
eled must be larger than the individual contribution of 
any of the stars included in the model. Using Figure 2 
from Cervino & Luridiana (2004) for Z= 0.4Z©, we com- 
pare the absolute magnitudes of the associations with the 
LLL in both the V and I bands. For young star clusters 
(~ 5 Myrs), the minimum V and I band luminosities for 
which population synthesis models are valid are -10.0 and 
-10.5, respectively. In the UV, the LLL would be fainter, 
and the effects of undersampling will be less. Since the 
associations have absolute V and I band magnitudes sig- 
nificantly fainter than the limits, it is unlikely we are 
completely sampling the stellar population. Stochastic 
effects then will impact the data interpretation, and we 
explore these effects more thoroughly in Section [3?3l 
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3.2. Independent Calculation of Cluster Ages and 
Masses 

We further constrain the stellar populations of the as- 
sociations by performing an independent estimate of the 
total number of ionizing stars present in each associa- 
tion. We extrapolate this calculation to determine total 
masses and ages of the associations. We begin by con- 
verting an Ha luminosity to an emission of ionizing pho- 
tons, Qo — 7.33 X 10^^ Lj^Q, s~^ with Lho in ergs s~^ 
cm-2 ^-1 ()Osterbrocklll989D . Models then predict the 
number of ionizing photons released per unit time from a 
certain spectral type (ST) of a massive star, usually 07V. 
By dividing the two quantities, one arrives at the number 
of "equivalent" stars of the given subtype. This number 
of "equivalent" ionizing stars, while a simple character- 
ization of the stellar population, does not represent the 
total number of ionizing stars. To determine the total 
number of ionizing stars, one must integrate over an ini- 
tial mass function, and use models that predict how the 
output of ionizing photons varies with stellar mass. In 
this section, we give the results of our calculation. In the 
appendix, we include a detailed, step-by-step description 
of the calculation, including all relevant equations and 
parameters. Ess entially , our m ethod follows the calcula- 
tions outhned in lVaccal ()1994f ). but uses t he recently up- 



dated Te// - spectral type calibration of iMartins et al.l 

dlpos). 

For associations 1, 2, and 5 we determine the upper 
limits of the present-day mass function to be 35 Mq, 25 
Mq, and 40 M©, respectively. The method we use, as 
described in the appendix, compares the total UV and 
Ha fluxes of the associations and the relative contribu- 
tions from O and B stars. Essentially, we are modeling 
the UV flux as the sum of the flux from both the O and 
B stars. We ultimately determine what value of stellar 
mass corresponds to a present-day mass function upper 
limit that can account for both the UV and Ha fluxes 
that are observed. We use these upper masses to cal- 
culate a limit on the lifetimes of the clusters. For initial 
stellar rotation velocities between and 300 km s~^, a 35 
Mq star burns hydrogen for 4.7 - 5.9 Myr, a 25 M© star 
burns hydrogen for 5.9 - 7. 4 Myr, and a 40 M© star b urns 
hydrogen for 4.2 - 5.1 Myr (jMevnet fc Maederl[200l . We 
obtain the total masses of the associations by integrat- 
ing the present-day mass function, and using the upper 
mass cut-offs given above and lower mass limit of 0.5 M©. 
The masses of associations 1,2, and 5 calculated by this 
method are 2700, 7200, and 970 M©, respectively. We 
note that this method of calculating the masses does not 
assume a scaled-down full IMF, as the Starburst99 mass 
estimate does, but rather uses a present-day mass func- 
tion. We integrate only to the present-day upper mass 
limit, resulting in a more physical calculation of mass, 
minimizing the number of fractional stars. The number 
of "equivalent" 07V stars, total number of ionizing O 
stars, the total mass, and the derived ages from both the 
Starburst99 model and this independent calculation are 
given in Table 5 for each isolated H II region. To ab- 
breviate the names of the two methods, we refer to the 
Starburst99 derived parameters with S99 , and t he ca l- 
culation using the i iiethod and models o f IVaccal (|1994D : 
IVacca et al.l (|l996D : IMartins et all (|2005f) as VM05. The 
range in ages given in Table 5 for VM05 correspond to 



stellar lifetimes without rotation, and the minimum error 
contributed by the uncertainty in the photometry. 

There is general qualitative agreement between the cor- 
responding quantities, yet some differences are clear. The 
ages determined with the VM05 calculation are system- 
atically larger than those determined from Starburst99. 
Neither model accounts for stochastic effects upon the 
stellar masses present, so the differences in ages and 
masses most likely reflect the difficulty of modelling the 
integrated properties of an undersampled stellar popu- 
lation. Association 2 is the most discrepant across the 
separate calculations, particularly in mass. 

A comparison between the number of equivalent 07V 
stars and the total O star population is illustrative of the 
limited information gained from a simple representative 
population. While association 1 has the most "equiv- 
alent" 07V stars, association 2 actually has a greater 
total O star population. This can be understood as fol- 
lows: association 1 is younger than association 2, and 
thus contains more massive O stars with stronger ioniz- 
ing fluxes than association 2; association 2 has evolved 
more than association 1, its most massive O stars having 
died out, and it is left with a large late- type O and B star 
UV radiation-producing population. This result is also 
readily apparent in the masses of the two associations. 

3.3. Stochastic Effects 

As with the interpretation of results from the Star- 
burst99 model, we caution that stochastic effects may 
play a large role, largely unaccounted for in the stellar 
population analysis of these associations. In the VM05 
calculation we have again assumed a fully populated 
IMF, which is most likely not the case, given the limits 
discus sed above and tabulated in ICervino &: Luridianal 
dlooi). StiU, we may probe the extent to which we are 
overestimating the stellar populations by examining the 
behavior of the Salpeter IMF we have assumed thus far. 
We assume that the IMF is a probability distribution 
function such that the mass of the largest star formed 
in an underpopulated cluster of N stars decreases as N 
decreases. The expression for this expectation value of 
the maximum mass star, < rrimax >, derived by Oey and 
Clarke (2005) is: 
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In this equation, $(m) is the IMF, integrated over 
mass, shown as two separate variables {m and M) for the 
two separate integrals. For a single power-law Salpeter 
IMF and M^p of 120 M©, the expected maximum mass 
star produced does not approach M^p until N stars with 
masses above 10 M© reaches ~ 1000 (in actuality, the ex- 
pectation value never reaches M^p as the integrated func- 
tion is asymptotic). Below this value of N stars, varying 
degrees of stochastic effects must be considered. Since 
we have already normalized the IMF for each associa- 
tion using the UV-flux from B stars, we perform another 
simple integration to calculate the number of stars with 
masses initially greater than 10 M© in each association. 
Association 1 contains 25 stars greater than 10 M©, as- 
sociation 2 contains 60, and association 5 contains only 
9. We numerically integrated equation 1 for N =25, 60, 
and 9, and for a Salpeter IMF to arrive at < nimax > 
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for each of the associations. These values are 80 M©, 
96 Mq, and 55 Mq for associations 1, 2, and 5. These 
results indicate that our age estimates and total masses 
for association 5 may be over estimated - if the maxi- 
mum mass star in a population of this size is initially 55 
Mq , and we have calculated the present-day upper mass 
limit to be 40 Mq as well, we could just as easily be wit- 
nessing a zero-age population forming in-situ. Stochas- 
tic effects are less dramatic for associations 1 and 2, but 
nonetheless important. The greater difference between 
the expected initial maximum mass and the present-day 
upper mass limit indicates directly that some evolution 
has occurred for associations 1 and 2. We stress that 
while our age estimates for association 5 are strict upper 
limits, we can have more confidence in our age estimates 
for associations 1 and 2. 

4. DISCUSSION 

In this section we discuss the properties of the associ- 
ations, clumps, and blue field objects, and compare and 
contrast them with those of stellar populations in the 
Galaxy and intergalactic H II regions of other systems. 
Do stellar clusters carry an imprint of the environment 
in which they form, or are their properties independent 
of their surroundings? We also examine the H I proper- 
ties of the environment containing the isolated associa- 
tions and comment on the conditions for their formation 
and potential evolution. Finally, the age and mass con- 
straints, along with physical sizes and locations, provide 
additional clues as to their likely fates in light of the rapid 
cluster dissolution inferred in many systems. 

4.1. Comparison With Galactic OB Associations, Open 
Clusters, and Other Intergalactic H II Regions 

4.1.1. Associations 

The associations span roughly 100 pc, and contain 
components separated by up to 10 pc, similar to the un- 
bound OB associations of the Milky Way. In addition, 
their luminosities and masses match those of a sample 
of Milky Way associations within 1 kpc, and the Galac- 
tic phenomenon of similarly-sized association su bgroups 
is dir ectly comparable to the clumps seen here (jBlaauwl 
Il964f ). The properties of the associations match those 
of the extragalactic OB associations in seven nearby spi- 
ral galax ies with dist ances ranging from 6.5 - 14.5 Mpc 
(jBresolin et al.lllQQSl ) equally as well. These extragalac- 
tic OB associations (diameters < 200 pc) have median 
diameters between 70 - 100 pc, contain 6 -15 blue stars 
with My < —4.76, and contain cluster candidates with 
V-I colors between -0.07 and 1.53. By comparison, the 
associations discussed here contain 3-9 clumps with 
My < -3.4, and V-I colors between -0.90 and 0.2. Pre- 
vious studies of spiral galaxy OB associations have noted 
a constant mean size of OB associations near 80 pc, a size 
that may represent the smallest step in t he hierarchical 
structure of giant star-forming regions (|Efremov et al.l 
[19871 : lEfremovl[l995l :l! vano ^ [l996h . The sizes of these in- 
tergalactic OB associations are consistent with the mean 
sizes of both Milky Way and extragalactic OB associa- 
tions. 

As noted in the Introduction, there is a growing num- 
ber of objects in the literature termed "intergalactic 
H II regions," with photometric properties that span 



a wide range, and environments that range from com- 
pact groups of galaxies to tidal streams from interact- 
ing galaxy pairs. One population of four intergalactic 
H II regions lies in the H I tail to the east of NGC 
7319, a member of Stephan's Quintet, a Hickson Com- 
pact Group composed of 3 intera cting galaxies at 80 Mpc 
^Mendes de Oliveira et al.ll2004f) . The four intergalactic 
H II regions in Stephan's Quintet appear to represent 
brighter, more massive, and larger versions of the inter- 
galactic H II regions of NGC 1533. Their absolute B 
band magnitudes cluster near -12, at the level of a faint 
dwarf galaxy, their derived masses are near 3x 10"* M©, 
their physical sizes (FWHMs) are between 400 - 600 pc, 
and their Ha luminosities are between 10^^ and 10'^^ ergs 
s^^. In contrast, the Ha luminosities of the NGC 1533 
intergalactic H II regions arc between 10'^^ and 10"^^ ergs 
s~^, with physical sizes < 100 pc, and masses on the 
order of 10^ Mq. 

In the ring-like H I structure surrounding the galaxy 
NGC 5291 (d = 62 Mpc), 29 intergalactic H II re- 
gions have been discovered, and appear to be charac- 
terized by a "UV-exc ess" compared to their Ha emission 
(jBoquien et al.l[2007t l. The ratio of the Ha star forma- 
tion rate (SFR) to the NUV SFR (SFR(Ha)/SFR(NUV) 
~ 0.5 for their entire sample) is considerably lower than 
typical ratios for Galactic H II regions (1 - 3). Here, 
we calculate SFR(Ha)/SFR(NUV) for our three inter- 
galactic H II regions, transforming our F250W UV fluxes 
(2715A) to NUV fluxes (227lA) by assuming a UV spec- 
tral slope typical of starbursting systems, F^ ^ A^^'^, 
thus multipl ying the UV fluxes given in Table 2 by a fac- 
tor of 1.56 (TMeurer et al.lfl9 991. We use the equations 
presented in |Bpquien et al.' (2007) for the SFRs (con- 
verted to cgs): SFR(UV) = [NUV] = 2.41 x lO"*" Lnuv 
[ergs s-^A-i] M0 yr"!; and SFR(Ha) = [Ha] = 7.9 x 
10~^^ Lho [ergs s"-"^] Mq yr^^. Associations 1, 2 and 5 
have [Ha] /[NUV] of 3.0, 1.3, and 3.1, respectively, and 
thus do not show any UV excess as observ ed for the inter- 
galact ic H II regions around NGC 5291. iBoquien et al.l 
(|2007t l explain the UV-excess in their sample as the re- 
sult of a quasi-instantaneous burst of star formation, re- 
sulting in two populations of stars, one roughly 4 Myrs 
old, and the other 10 Myrs. For comparison, their PSF 
in the NUV has a characteristic width of 5", 75 times 
larger than our PSF. At 62 Mpc, 5" covers 1.5 kpc, and 
would therefore encompass most of the HRC field of view 
at the distance to NGC 1533, and quite possibly all of 
the visible star formation around regions 1 and 2, which 
we show spans a range of possible ages (2-9 Myrs). 

4.1.2. Clumps and Blue Field Objects 

Aside from single O and B stars, the closest Galac- 
tic analogs to the clumps and blue field objects are 
young open clusters. For example. iBonatto et al.l ()2006D 
presents a detailed analysis of the properties of NGC 
6611 (the Eagle Nebula) deriving an age of 1.3±0.3 Myr, 
a limiting radius of 6.5±0.5 pc, and an observed mass of 
1600 Mq. Viewed from a distance of 19 Mpc, a young 
open cluster such as NGC 6611 would appear unresolved 
in the HRC and similar to one of the clumps or blue field 
objects. Moreover, the integrated photometric charac- 
teristics of Galactic open star clusters younger than 10 
Myrs are fairly consistent with those of the clumps and 
field objects. Of the 140 open clusters (masses between 
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TABLE 5 

Derived Properties of the Stellar Associations^ 



S99: VM05: 
HII No7V Mtot Age Nov ^tot Age 

Region 10^ Mq Myrs 10^ M© Myrs 



1 8.06 1.24-1.49 2.9 - 3.5 11 2.73 4.4 - 5.0 

2 5.46 1.85-2.03 4.2 - 4.8 22 7.17 5.6 - 6.2 
5 3.81 0.61-0.84 2.2 - 3.8 4 0.97 3.8 - 4.6 



^ No7V is the number of equivalent 07V stars in each association, calcu- 
lated using the SINGG Ha luminosity; Nov is the total number of O stars 
present in each association (see Section 3.2 and the appendix for details); 
Total stellar masses and population ages are given for two analysis methods, 
a comparison with Starburst99 model outputs (S99), and an independent 
calculation described in Section 3.2 (VM05). 



25 and 4 x 10* Mq) tabulated in lLata et all (|2002f ). 17 
are younger than 10 Myrs, with absolute V band magni- 
tudes between -4.5 and -9, and V-I colors (for only 6 of 
the 17) between and -0.27. These colors are similar to 
those found for the clumps and blue field objects listed 
in Tables 2, 3, and 4. However, the properties of the 
clumps and blue field objects are also in line with single 
O and B stars, seen in Figures 5 and 6. In the end, we 
cannot distinguish between single stars and young open 
clusters at the physical resolution limit of the HRC. 

The clumps and blue field objects also have extragalac- 
tic analogs, and more extreme counterparts. The photo- 
metric properties and sizes of the clumps and blue field 
objects align with those of the faint 1- 10 M yr c luster s 
in the LMC and SMC studied bv [Hunter et all (|2003l ). 
Compared to the young clusters in the disk of M82, a 
violently starbursting system, the clumps and blue field 
objects have similar V-I colors (on the blue side) and ab- 
solute V-band magnitudes. The range of V-I colors for 
the M82 young clusters is -0.05 to 2.34, with absolute 
y-band magnit udes ranging from -11.2 to -5 .3, having a 
mean of -6.75 (jKonstantopoulos et al.l [20071 ). However, 
the young massive 10^ M©) clusters of the Antennae 
galaxies exhibit a range of properties nearly on a differ- 
ent scale from those of the stellar populations powering 
NGC 1533's intergalactic H II regions. The clumps lie 
on the faintest end of the Antennae cluster lu minosity 
function constructed bv lWhitmore et aTl (|1999f ). contain 
far fewer stars (indeed, they may be single stars), and for 
the most part, occupy the bluest end of the V-I color dis- 
tribution. Overall, the properties of the clumps and blue 
field objects show no deviation from the continuum of 
properties for known star clusters (if they are not single 
stars) in both our Galaxy and others. 

4.2. Gas Properties 

The H I distribution shown in Figure 1 reveals a ring 
around NGC 1533 with a total H I mass of 6 x 10^ Mg 
(jRvan- Weber et all [2003bl ). The southeast cloud, the 
site of the star formation discussed throughout this pa- 
per, contains approximately 1/3 of this total H I mass 
(2.4 X 10^ Mq), and has a maximum H I column den- 
sity of 4 X 10^° cm~^. The three associations do not lie 
in the densest regions of the SE H I cloud, but rather 
in regions well below the canonical neutral gas surface 
density cutoff for star formation on global scales (10^^ 
cm^^, Skillman 1987). Of course, the Q scenario pre- 
dicts that this critical gas density threshold (Scrit) scales 



with velocity and radius (jKennicu ttl ll989f). For a flat ro- 
tation curve and parameters fitted in iKennicuttj (|1989D . 
^crit (Mq pc~^) = 0.4 where c is velocity disper- 
sion and V is the rotation velocity of the galaxy, both in 
km s^^, and R is the projected distance from the galaxy 
nucleus in kpc. Given a fitted rotation velocity of the 
gas around NGC 1533 of 150 km s~^ and a velocity dis- 
persion of 10-30 km s-i (|Rvan- Weber et al.|[2003a[ ). at 
the projected radii of the isolated H II regions (19 - 31 
kpc), the threshold column densities of the H I would 
range between 4 — 20 x 10^° cm~^. The observed col- 
umn densities at the positions of the intergalactic H II 
regions appear to fall below this range in critical surface 
density (see Table 1); although it is possible the true H I 
column densities are underestimated due to the 1 arcmin 
beam. It should also be noted that though the H I ring 
appears to be rotating around NGC 15 33, it is not part o f 
a smoothly rotating disk as studied bv lKennicuttI (|1989f) . 

Across entire galaxies, star formation regions tend to 
align with H I peaks, but this g lobal correlation breaks 
down on sub-kiloparsec scales (jKennicutt et al.l[2007| ). 
On a local scale, low-levels of star formation can pro- 
ceed despite lower than normal neutral gas column den- 
sities. Similar to the H II regions outside NGC 1533, 
the Stephan's Quintet intergalactic H II regions are not 
associated with H I peaks in the tidal tail (at le vels be- 
low 10^° cm-2) (jMendes de Oliveira et ani2004l ). Addi- 
tional ly, the H II region found in the Virgo Cluster by 
Gerh ard et all (2002 ) lies in no detectable amount of H I 
gas (down to lO^^cm"^; Oosterloo & Van Gorkom 2005). 
There must be a reservoir of yet undetected molecular 
material at the positions of these intergalactic H II re- 
gions. The Atacama Large Millimeter Array (ALMA) 
will shed light on the molecular conditions required to 
trigger star formation outside the disk of galaxies. 

Though it may appear as though associations 1 and 2 
lie near a relatively dense cloud of H I (just to the north 
of the associations), there are actually multiple clouds 
along the line of sight, and not a single, rotating cloud. 
Figure 8 shows the lack of a smooth velocity gradient 
at the position of this apparent H I concentration in a 
position velocity slice across the long diameter of the 
cloud. The total H I mass of the concentration of gas 
clouds in the vicinity of associations 1 and 2 is 3 x 10* 
Mq, similar to the amount of H I found in the vicin- 
ity of the potential t idal dwarf galaxie s of the Antennae 
(iHibbard et al.l[200ll) and NGC 5291 (|Bournaud fc Dud 
l2006f l. but lower than what is typically found for tidal 
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Fig. 8. — A slice in position- velocity space across the long diam- 
eter of the apparent H I concentration just to the north of associa- 
tions 1 and 2. This figure shows that the concentration is made up 
of numerous clouds across a wide range of velocities. The noise level 
isat~ 4 mjy/bm and the beam is 1 arcminute tRvan-Weber et al.l 
I2003bl ). 

dwarf galaxies (> 10'^ Mq e.g. Due & Mirabel 1994, 
1998; Hibbard & Barnes 2004; ). The velocity disper- 
sions in the vicinity of NGC 1533's intergalactic H II 
regions are near 30 km s~^, with gra dients in the range 
7-50 km s-i kpc-i ()Rvan- Weber et al. 2003a). The rel- 
atively high velocity dispersions may lend some clue to 
the conditions of their formation, as cloud-cloud colli- 
sions are known to trigger star formation. These velocity 
dispersions are much higher than the H 1 velocity disper- 
sions measured in the vicinity of the tida l dwarf galaxy 
candid ates in the Antennae (5-15 km s~^: iHibbard et all 
()2001h ). and the overall kinematic structure is much more 
irregular than that of both the Antennae and other po- 
tential tidal dwarf galaxies (e.g., Bournaud et al. 2007). 
We plan to explore the relation between H I column den- 
sities and velocity dispersions for star formation in the far 
outskirts of galaxies in a future paper including a larger 
sample of intergalactic H II regions and high-resolution 
H I data. 

While the stellar masses revealed here do not place the 
star formation in the HI ring in the category of tidal 
dwarf galaxy (TDG) candidates, the location of associ- 
ations 1 and 2 in the tip of the SE H I tidal tail may 
actually be a pri me one for surviva l . In a set of 96 N- 
body simulations. iBournaud fc Dud ()2006f ) found that of 
all structures formed during mergers, the only objects 
that are long-lived are those formed in the extremities of 
the tidal tails. The tip s of tails are ofte n sites of kinemat- 
ical pilings- up of gas (' Due et ani2004f ). contain popula- 
tions less prone to fall back upon the progenitor galaxy, 
and lie far beyond the r each of disruptive galactic tidal 
forces (jBournaud fc Dug .2006.) . However, these simula- 
tions were limited to stellar tidal structures with masses 
greater than 10* Mq, a criterion our associations surely 
do not meet. If one assumes the gas in the clump near 
associations 1 and 2 is part of a bound structure, the 
resulting dynamical mass is 20-25 times higher than the 
baryonic mass in this region. This abnormally high dy- 



namical mass for a TDG, together with the structure of 
the clouds, indicate the stellar associations are unlikely 
to evolve into a relatively long-lived tidal dwarf galaxy. 

4.3. Fate and Evolution 

In the recent literature, there has been much focus on 
the age distribution of star clusters as a way of exam- 
ining their rate and nature of dissolution. We cannot 
construct an age distribution of the 19 individual clumps 
and 18 field objects because they all have F-band ab- 
solute magnitudes considerably fainter than -9, a limit 
indicating that single stars dominate their output radia- 
tion. We can, however, comment on the associations and 
their relationship to the clumps and blue field objects in 
the context of cluster disassociation. In this section, we 
assess the survival of the associations by considering the 
locations of the objects in the HRC fields relative to each 
other, their relative sizes, and age and mass constraints 
from stellar evolution models. 

4.3.1. Rapid Dissolution 

Star clusters of a wide range of masses and in vari- 
ous environments are subject to a large number of pos- 
sible internal and external evolutionary effects. In the 
Milky Way, embedded protoclusters less massive and 
more densely concentrated than the isolated associations 
discussed here are observed to dissolve on fast timescales 
("infant mortal ity"), with less tha n 4-7% surviving as 
bound clusters ()Lada fc Ladall20d3l ). In the extreme en- 
vironment of the colliding Antennae galaxies, the sur- 
vival prospect is just as grim for mo re massive young 
"super star" clusters ([Fall et al.ll2005D . In more quies- 
cent systems like the SMC, "infant mortality" may play 
a much smaller role, if any; however, its role has not 
been fully investigated for clusters younger than ~10 
Myrs due to a detection bias against c lusters having neb- 
ular emission in the SMC samples (jCieles et al.l l2007t 
Ide Grijs fc GoodwinI l2007t ) . The fate of the stellar sys- 
tems powering intergalactic H II regions (ages < 10 Myrs; 
masses 10'^ - lO'* M©) has not yet been explored in the lit- 
erature. On one hand, their fate may be similar to that 
of Galactic OB associations, thought to expan d freely 
for tens of Myrs (|Blaauwlll96l Ide Bruiindll9990 . In this 
case, one of their most massive clumps (subgroups) on 
average e nds up as a bound, open cluster such as the 
Pleiades fcada fc Lada 2003). On the other hand, their 
location in an environment vastly different from that of 
the Galactic disk may result in a different fate. Further- 
more, objects of similar masses and sizes to these isolated 
associations seem to be able to s urvive for hundreds of 
Myrs in the SMC ()de Griis fc Go odwin 2007), while the 
vast majority appear to dissolve o n the order of tens 
of Myrs in the Antennae Galaxies ((Fall et al.ll2005f) . In 
this context, it is relevant to explore what will become 
of these isolated associations given our high-resolution 
photometric data. 

The scenario of rapid cluster dispersion requires some 
disruptive effect other than the general stellar evolution- 
ary processes and gravitational effects such as stellar 
evaporation driven by two-body relaxation considered in 
typical cluster fading laws. Most commonly proposed is 
that the energy and momentum input from massive stars 
(in the form of winds, jets, ionizing radiation, and super- 
novae) remove the cluster's ISM, leaving it unbound and 
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freely expanding piilifToSOl : lGoodwinlfT997t iFall et all 
l2005f ). These effects result in an average expansion rate 
of the stars in Galactic OB associations o f 5 km s~^ , 
or 5 pc in 1 Myr, estimated analytically (jHillsl [l980l ) . 
Obseryationally, the direct measurement of an expan- 
sion speed will haye to wait until the launch of Gaia. 
Howeyer, under realistic assumptions of expansion, the 
stars' proper motions and radial velocities (along with 
the best positions available from Hipparcos) are con- 
sistent with linear expansion coefhcients of about 5 km 
s-i per 100 pc tBlaauw..l964 : .Bobvlev fc B aikova 2007il. 
Thus, typical Galactic associations can expand for only 
tens of Myrs before dispersing into the general field pop- 
ulation. If these intergalactic OB associations expand 
similarly, they eventually will not be recognizable as dis- 
tinct groups of stars. The expansion of associations also 
implies that larger associations, will be, on average, older 
than more compact associations. Of the three studied 
here. Association 5 is the youngest (2 - 4.5 Myrs), and 
also has the smallest visible physical extent of roughly 
60 pc. On the other hand, association 2 is the oldest 
(4.5 - 6.5 Myrs), and has the largest physical extent of 
roughly 100 pc. This age-size sequence is consistent with 
a scenario of dispersion. Association 5 may still be recog- 
nizable as an association of massive stars for some time, 
whereas associations 1 and 2 likely have less time before 
they will be indistinguishable from the blue field objects. 

Do the blue field objects then represent stars and/or 
star clusters that were once physically linked with the as- 
sociations and have since dispersed into the field? Given 
the typical expansion rate of OB associations, if the blue 
field objects were once part of association 1 or 2, they 
would have to be 30 - 240 Myrs old to account for the 
observed separations of 150 - 1200 pc. The colors of the 
blue field objects are not consistent with this range of 
ages. The single stars dominating the radiation output 
of the blue field objects are roughly consistent with main 
sequence stars ranging from 8 Mq to 40 Mq (consistent 
with B2V - 05V stars; see Figures 5 and 6). Geneva 
stellar tracks (iLeieune fc Schaeredl200l1 ) predict that an 
8 Mq star has a V-I color below zero as it burns hydro- 
gen, and reddens by over 1 magnitude at an age of 27 
Myrs, while a 40 Mq star has a V-I color below zero as 
it burns hydrogen, and reddens by more than 1 magni- 
tude near 6 Myrs. Since the blue field objects generally 
have V-I colors below zero, they are unlikely to be older 
than 6-27 Myrs, and therefore would not have had 
time to reach their current projected separations from 
the nearest association (1 or 2) at the assumed velocity 
of 5 km s"""^. There is also no indication that the field 
objects further from the association are redder (and, by 
inference, older) than those that are nearby. Some of 
the blue field objects may also represent single runaway 
O or B stars ejected from the associat ions at maximum 
speeds ranging from 200 - 350 k m (jHoogerwerf et al.l 
120011: [Leonard fc DuncanI [l990l ) . This speed would give 
the stars the ability to travel between 1200 - 9450 pc in 
their estimated lifetimes, but this phenomenon is highly 
unlikely to account for all 18 blue field objects, if any at 
aU. 

A more likely scenario for the blue field objects in the 
context of cluster dissolution is that some of them once 
formed their own association that is no longer easily rec- 
ognizable as such. In the Galaxy, the Cassiopeia- Taurus 



group may represent a similar grouping of objects with a 
common origin in a more advariced st ate of disintegration 
(|Blaauwlll96l Ide Zeeuw et al.lll999[ ). Here, F5, F6, F8, 
F9, and Fll (and other tiny sources visible in the same 
locale, but undetected by SExtractor, see Figures 2 and 
4) appear to be clustered together themselves, in an os- 
tensibly young, diffusing complex with a total extent of 
about 150 pc and average object projected separations 
of 20 pc. The radius of this potential complex, given av- 
erage association expansion velocities and assuming all 
its components formed at roughly the same position, is 
indicative of an age of ^ 20 Myrs, which is roughly con- 
sistent with the range of ages possible for the blue field 
objects (see above). In the SINGG Ha image (1.37" res- 
olution), we do not detect unresolved Ha emission at the 
location of this "complex" as we do for associations 1, 2, 
and 5. Assuming the H II region potentially associated 
with this diffusing complex is similarly unresolved in the 
SINGG image (likely, given the 1" FWHM of the com- 
plex of objects), the upper limit to the Ha flux at the 
location of this "complex" is 9.96x10"^^ ergs cm~^, 
based on 5a point source detection hmits of the SINGG 
data. Dividing this upper limit by the sum of the UV 
flux of F5, F6, F8, F9 and Fll, we derive an upper limit 
Log {Fho/F^^^^^) of 1.44. The subsequent lower age 
limit of the diffusing "complex" based on Figure 7 is 5.6 
Myrs, at least 1-2 Myrs older than the ages of associa- 
tions 1, 2 and 5 based on this same plot, and consistent 
with the idea that these blue fleld objects represent an 
association in the process of dissolving while the other 
associations are comparably more intact. 

One can also consider what should be observable in 
the field if we assume a continuous and homogeneous 
low rate of star formation in this region for the past 
100 Myrs. In this scenario, if we see 2 associations of 
2000 Mq that are under 10 Myrs old, then it follows 
that we would expect 2 associations of 2000 M© every 
10 Myrs for the last 100 Myrs, or twenty 2000 Mg as- 
sociations in this region. Independent of the processes 
responsible for rapid cluster dissolution, we should con- 
sider the effects of fading due to mass loss from stellar 
evolution and tidal effects, as this will decrease the num- 
ber of clusters we expect to see. Because our d etection 
limits are similar to those of lHunter et al.l ()2003l ) in their 
study of the Large and Small Magellanic Clouds, we as- 
sume a cluster fading function similar to theirs, where 
Mfade = 982(t/Gyr)°-69 Mq. Thus, without rapid clus- 
ter dissolution, we would be able to detect similar mass 
associations with ages up to ~1 Gyr (approximately the 
age of the HI ring, Ryan- Weber et al. 2003a ). That 
we do not detect associations independent of those di- 
rectly associated with the H II regions indicates that the 
clusters are either rapidly dispersing or this is the first 
generation of star formation in the HI stream. 

In the scenario of rapid dissolution, we may still expect 
to find background fluctuations in our images represen- 
tative of previous star formation in the HI stream that 
has subsequently dispersed. In our HRC data, we see no 
evidence of an underlying older population after investi- 
gating the noise distributions in each band. The shape 
of the noise distribution is nearly perfectly Gaussian in 
all cases, and in UV, V, and I we probe down to 28.0, 
28.9, and 29.0 AB magnitudes, respectively. While this 
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does not rule out an older population of dispersed clus- 
ters (or a previous stellar component to the HI stream), 
if one exists, it is below the noise level. Considering 
an expansion rate of 5 km s~^, we can estimate the ex- 
pected cluster fading due to its dissociation alone. A lin- 
early increasing cluster radius with time (|Blaauwlll964l : 
iFall et an2005f) . results in the surface brightness decreas- 
ing with the square of time (in fact, the surface bright- 
ness will decrease faster because as stars age and die, 
the cluster becomes fainter) . For areas the size of these 
associations (~186 pixels^), our HRC images give a 5a 
surface brightness limit of 24.3 mag arcsec"^ in the V 
band. We calculate the y-band surface brightnesses of 
the associations 1, 2, and 5 (using their half-light radii) 
to be 21.85, 22.26, and 22.0 mag arcsec"^, respectively. 
Thus, in 14 Myrs, when the average radius of the asso- 
ciations has increased by a factor of 2.5 (from 47 pc to 
117.5 pc), the surface brightness of the associations will 
fade below the detection limit of the HRC images. It is 
consistent with the idea of rapid cluster dissolution that 
we see one nearly dissolved association roughly 6-20 
Myrs old, and several younger, more intact associations. 

4.3.2. Alternative Scenario: A First Generation of Stars 

Contrary to the above argument for rapid dissolution, 
the associations, clumps, and blue field objects may rep- 
resent a distribution of a first-generation of star forma- 
tion in which the stars and clusters have not strayed far 
from their location of formation. We emphasize that 
this scenario is self-consistent only if star formation in 
this region began less than 10 - 15 Myrs ago, otherwise 
we would see older, intact associations similar in mass 
to those observed. In this scenario, the blue field objects 
would occupy the low mass end of the initial cluster mass 
function (CMF) rather than represent the remnants of 
dissolved larger stellar associations. That the youngest 
association (association 5) is the most compact may sim- 
ply indicate a lower initial mass of this association com- 
pared to the others, and not that it is part of an age-size 
sequence of cluster dissolution. If the observed associa- 
tions comprise a first generation of stars, then we cannot 
use these HRC data to untangle their fates. 

In this regard, the isolated and singular association 5 
may be most likely to be an initial episode of star forma- 
tion. For, if star formation has been continuous in this 
region for any period of time greater than 10 Myrs, the 
absence of any sign of dispersed cluster remnants in its 
vicinity may be inconsistent with cluster disintegration. 
That this youngest association remains isolated in com- 
parison with its older counterparts offers some indication 
that it represents an initial episode of star formation in 
this northern region of NGC1533's H I ring. There are a 
few related caveats to consider: 1) even if association 5 
represents a first generation of star formation, rapid dis- 
solution may subsequently proceed 2) association 5 could 
simply represent the youngest, most compact component 
of the age-size sequence we expect for cluster dissolution 
3) the star formation in this northern region may have 
begun later than that in the southeast, and so should 
simply be considered separately. 

5. SUMMARY AND CONCLUSION 

HST ACS/HRC images in UV, V and I bands resolve 
the three very young stellar associations powering the 



intergalactic H II regions that lie up to 30 kpc from the 
stellar disk of the SO galaxy NGC 1533. We have cat- 
alogued and performed circular aperture photometry on 
48 individual sources in the HRC images, 27 of which 
appear to be associated with the ongoing star formation. 
From background galaxy statistics and broadband colors, 
we have determined that the other 21 sources are likely 
distant background galaxies. Furthermore, we have de- 
fined a range of sources in the images: the term "asso- 
ciations" refers to the total stellar populations powering 
the H II regions; the term "clumps" refers to the indi- 
vidual resolved sources composing the associations; and 
the term "field objects" refers to the individual sources 
detected separately from the associations. Using two 
different methods, (Starburst99; Leitherer et al. 1999, 
and an independent calculation derived from Vacca 1994 
and Martins 2005, see appendix) in combination with our 
broadband co lors and magnitude s, and Ha luminosities 
from SINGG (jMeurer et al.ll200^ . we have placed upper 
limits on the the total numbers of stars, masses and ages 
of each association. The masses of the isolated associa- 
tions range from 600 Mq to 7000 Mq , and the ages range 
from 2-6 Myrs. Given the faint luminosities of the stel- 
lar associations, we have also considered the importance 
of stochastic effects upon our calculations. Colors and 
absolute magnitudes of many of the clumps within the 
associations and field objects are consistent with single 
O and B stars dominating the output radiation. 

There are three results from our analysis that deserve 
emphasis: 

• Despite their location in an environment that is 
vastly different from that of a Galactic spiral arm, 
the photometric properties and sizes of the associa- 
tions, clumps, and blue field objects are consistent 
with those of Galactic stellar associations and clus- 
ters and/or stars. This result suggests something 
uniform about the star-formation process on scales 
of tens of parsecs, independent of environment. 

• The level of star formation in the SE H I cloud 
of NGC 1533 is lower than that observed in other 
systems (e.g. NGC 5291; Boquien et al. 2007). 
Furthermore, low H I column densities (lO^^cm^ ) 
and high H I velocity dispersions (~ 30km s^^) 
characterize its environment. These data hint that 
there may be some relation between SFR, H I col- 
umn density and velocity dispersion that holds out- 
side of massive gas disks, on local scales. These 
data are also consistent with the recent finding 
that H I surface density and star formation rate 
densities do not correl ate on sub-kiloparsec scales 
(jKennicutt et al.l[2007h . The kinematics of the H I 
in the vicinity of the H II regions do not suggest 
the stellar associations will evolve into a larger tidal 
dwarf galaxy. 

• The stellar associations detected here are either in 
the process of rapid dissolution or represent the 
first generation of stars in the gaseous ring around 
NGC 1533. If the star formation in the SE part of 
the H I ring has been continuous over the last 20 
Myrs to 1 Gyr, we are directly witnessing the rapid 
dissolution of young clusters and associations. Not 
only do the associations display a range of ages and 
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sizes consistent with cluster dissolution, but the 
age constraints and locations of several clustered 
blue field objects indicate they are the last visible 
remnants of an already dispersed association. If 
instead, these isolated associations and blue field 
objects represent a first generation of stars in the 
stripped H I of NGC 1533, we cannot determine 
the ultimate fate of the associations. This possi- 
bility of a first generation of stars is then uniquely 
suited for exploring star formation triggers and ini- 
tial conditions without the influence of stellar feed- 
back. Given that their Galactic OB association 
analogs are generally thought t o be gray itationally 
unbound and freely expanding ()Blaauw| [l964: Hills 
[198d), along with evidence for star cluster "infant 
mortality" in various clustered stellar systems (e.g. 
Lada & Lada 2003; Fall et al. 2005), we find it 
most likely that the stellar populations of the in- 
tergalactic H II regions will soon (< 20 Myrs or 
so) no longer be detectable nor recognizable as as- 



sociations of stars; rather, they will be building a 
background population of young objects in the halo 
of NGC 1533. 
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APPENDIX 

INDEPENDENT CALCULATION OF CLUSTER AGES AND MASSES 

Here we pr esent the det ails of the estimate of the total number of ionizing stars present in each star-forming region 
derived from lVaccal ()1994[ ). In Section 3.2, we use this method to determine association ages and masses, and compare 
the results with those of Starburst99. This method assumes no stellar evolution off the zero age main sequence, but is 
still provides reasonable estimates for the total number of ma in sequence O stars because the emission rate of ionizing 
photons does not vary substantially as massive stars evolve (|Vaccal[l994[ ). In A.l, we outline the calculations of the 
total number of O stars based on the observed ionizing flux of each association. In A. 2, we use the observed UV fluxes 
to constrain the number of B stars, set the upper mass limit to the present-day mass function, and normalize the IMF 
for each isolated association. 

The Number of Stars 

An Ha lumin osity corresponds to an emission rate of ionizing photons, Qo = 7.33 x 10^^ L//q s^^ with Lho in ergs 
cm~^ (|Osterbrocklll989l ). Using this conversion factor, one can derive the number of "equivalent" stars of a 
given subtype powering an H II region. This number of "equivalent" ionizing stars, while a simple char acterization o f 
the stellar population, does not represent the total number of ionizing stars. With this notion in mind, IVaccal ()1994[ ) 
present values for the conversion factor (770 = N^yy/Noy) between the "equivalent" number of O stars of a given 
subtype and the full ionizing population, utilizing the most current stellar atmosphere and stellar evolution models 
available in 1994. In this appendix, we describe our calculation of the total ionizing population of each H II region. 
Essentially, our me thod follows the calcu l ations o utline d in .Vaccal ()1994D . but uses the recently updated Teff - spectral 
type calibration of iMartins et al.l (|2005D . IVaccal (|1994i ) thoroughly details the calculation of this conversion factor, 770, 
and we refer the reader to that paper for a full description of the calculation. 

The motivation for recalibrating 770 using Martins et al. (2005) arose from the significant difference in their stellar 
parameters from previous results. Taking non-LTE, and wind and line-blanketing effects into account leads to effective 
temperature scales that are cooler than previously repor ted scales by 2000 to 8000K. Furthermore, the inclusion of 
metals in the model atmospheres of IMartins et al.l (|2005[ ) modifies the ionizing fluxes of the O stars. A true detailed 
recalculation of rjQ would have to include the new evolutionary tracks and regenerate stellar models, a task beyond the 
scope of this paper. Instead, we present a "re-estimation" of the conversion factor r/o that hinges upon the number 
of ionizing photons produced as a function of stellar mass, Qo(M). Figu re M shows this fun ctio n plotted versu s stellar 
mass for two sets of stellar parameters taken from tables pr esented in iVacca et al.l (|1996f) and IMartins et ah! (|2005i ). 
Below 15.8 Mq, we used ionizing fluxes read from figure 1 in IVaccal (|1994l ) for both plots, making our two power-law 
fits identical in this regime. This approximation is valid since 1) the lower mass stars contribute very little to the 
ionizing radiation and 2) the two T^ff - spectral type calibrations are roughly the same below about 16 Mq. The solid 
lines represent the broken power-law fits (a series of linear fits in log-log space) to the data. 

We begin our re-estimation of ? 7„ by considering the equation for the conversion f actor presente d in IVaccal (|1994f ) and 
applying the values tabulated in IMartins et al.l (|2005f ). From physical definitions. IVaccal (|1994( ) derives the following 
equation: 



Vo 



= N'orv/Nov 



Here, $(Af) is the present day mass function with the form AM~(^+^\ where A is the normalization constant and x 
is the power law index (x=1.35 for a Salpeter mass function). Q^"^^^'^(M) is the Lynian-contimmm emission rate in 
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Fig. 9. — On the left: The ratio of the ionizing photon emission rates in photons as a function of stellar mass between two sets of 
stellar parameters, those presented in Vacca et al. (1996) and Martins et al. (2005). The separate functions Qo(M) are quite similar, yet the 
differences produce significant changes in the conversion factors used to determine the total ionizing stellar populations (see Table 4 and the 
appendix for values and a detailed descri ption). Below 16 M©, the two functions are exactly th e same. On t he right: The polynomial fits 
to the data points of lMartins et all (|2005l ) are plotted. The masses used to ge nerate Qo(M) from I Vacca et aP ((19960 were the evolutionary 
masses, while the masses used to generate Qo(M) from [Martins et al.l ( 120051) were the spectroscopic masses (the only masses tabulated). 
Using the spectroscopic masses does introduce some error to our calcula t ion, t hough slightly tempered by our use of only those masses 
below 60 Mq, and the lower effective temperature scales of lMartins et all (120051) . 



0.070 




phot ons s ^ as a func tion of stellar mass for stars on the main sequence. Using the relatio n Qp — 7.33 x IO^^ Lho 
s-i (|Osterbrocklll989l ). the Ha luminosities given in Table 1, and the fiducial Q^'^^ = 48.63 (jMartins et al.ll2005l ). we 
calculate the number of equivalent 07V stars ( ^'o7v) associations 1, 2 and 5 to be 8.06, 5.46, and 3.81, respectively. 
We then evaluate Equation Al to calculate 770- From the definition 770 = NQ^y/Nov, we calculate the total ionizing 
main sequence O star population (Nov) fo r a range of upper in ass limits (Mup) given Qf"^*^'^(M) from Figure [9)3. 
The lower mass (M;o„,) is taken to b e 3 M0 ()Schmidt-Kaleijll982[ ) . and the lower mass limit to O stars {Mqv) is taken 
to be 15.8 M0 ([Martins et al.ll2"005f ). The total numbers of main sequence O stars in each association are estimated 
at 11, 22, and 4 for associations 1, 2 and 5 using the present-day mass function upper mass limits we determine in the 
following section. 

Table 6 gives three different calculations of ?7n for sola r metallicity, a Salpeter IMF, and various values of Mup. 
The first column of rjo is taken directly from 
comparison, we also give ?7o for the Qo(M) of 



Vacca 



Vacca 



(Il994l) . in which rjo was calculated at larger mass intervals. For 
([1994) calculated by our method of function fitting, integrating, 
and our choice of a lower mass for O st ars of 15.8 M ^. t^q is this re-estimation of Vacc a's original c onversion factor 
using the recent stellar models of Martins ' et al.l (j2005[ ) . Our values of rjQ for the model of lVaccar(|1994[ ) agree well with 
the values of tjq gi ven directly in Vacca (1994'). Furthermore, it is readily apparent that the most recent Te//-spectral 
type calibration of lMartins et al., (2005) has impacted the values of r/o quite significantly. 

We caution that these values of rjo should not be used as more than estimates, because of the ass umptions and 
appro ximations used in the calculations. Still, these are the most up-to-date estimates of their kind. iMartins et al.l 
()2005[ ) quote errors in their masses between 35 and 50%, which results in rjo having 15-20% errors. While this 
calculation technically requires the use of an evolutionary (model dependent) mass, Martins et al. give only the 

theoretical spectroscopic {Mspec = ) masses for each O star spectral type. Finally, the lower mass limit to O 
stars is quite poorly defined, which adds at least 10% to the error in the calculation of rjQ. The discrepancy between 
evolutionary and spectroscopic mas ses is a known i ssue w ith stellar models, but is somewhat alleviated by the decrease 
in the effective temperature scale of lMartins et al.l (|2005t ). We also stress that we have done this analysis for only main 
sequence solar-metallicity O stars. Varying the metallicities of the modeled stars will also have a small effect on the 
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outcome. 



The Number of B Stars 

We calculate the total non-ionizing B star population expected from the UV luminosity of each association taking 
a path similar to that taken above. If the total UV luminosity is produced by only O and B spectral types (to a good 
approximation), and we can quantify the ionizing population by the Ha luminosity, then we can also quantify the 
population of stars that contributes exclusively to the UV luminosity. Given the number of equivalent 07V stars from 
the Ha luminosity, we calculate the number of equivalent B9V stars from the UV luminosity as follows: 

Fuv = No7v X F^l^ + Nbw x (A2) 

In this equation, F^^^ is Fa of an 07V star in the F250W bandpass calculated using SYNPHOT (6.64 x 10"^^ ergs 
s-i cm-2 k-^ at the distance of NGC 1533 from Kurucz 1993 spectral atlas in the SYNPHOT database), F^^"^ is Fa 
of an B9V star in the F250W bandpass (9.20 x IQ-^z ergs s'^ cm-2 at the distance to NGC 1533), and Fuv is 
the total Fx of the object measured in the HRC F250W image. For associations 1, 2, and 5 we calculate 3000, 8200, 
and 1100 equivalent B9V stars respectively. 

Given the number of equivalent B9V stars for each association, an analysis similar to what we described above was 
undertaken to determine the total number of B stars that contribute to the UV luminosity of the intergalactic H II 
regions. In the same sense as for O stars, a conversion factor between a number of "equivalent" B stars of a certain 
spectral type and the full population of B stars can be determined: r]o{B) =N^gy / ^bv ■ In this case, N'^gy is the 
number of equivalent B9V stars (see above), and ^bv is the total B star population. The corresponding integral for 
the calculation of the conversion factor is simply: 

f ^^^-^ '^{M^F^^^^^ (M) dM 

We generate F^'^^^'^ (M) for all B spectral types using the masses and temperatures of ISchmidt-KaleJ ()1982f ) and 
calculating the flux in the F250W bandpass of each spectral type using SYNPHOT and the Kurucz spectral type 
models. $(Af ) is the same present-day mass function used previously. Next, we calculate the total number of non- 
ionizing B stars by integrating Equation A2 for masses between 3 and 15.8 M©. We estimate a total B star population 
of 150, 400, and 50 in associations 1, 2, and 5, respectively. 

Finally, for comparison with our original calculation, we use the normalization of the present day mass function 
from this calculation, and integrate $(M) between 15.8 and various values of M„p to determine the total number 
of main sequence O stars. The total number of O stars calculated using the Ha luminosity should match the total 
number of O stars calculated using the UV flux at some upper mass limit, interpreted as the maximum mass of the 
PDMF. For isol ated a s sociat ions 1, 2, and 5 these upper masses are 35 M0, 25 M©, and 40 Mq. Using the values of 
770 presented in IVaccal ()1994[ ) gives slightly higher upper masses of 40 Mq, 28 Mq, and 50 Mq for associations 1, 2, 
and 5, respectively. 
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TABLE 6 
New Estimates of the 

CONVERSION factor rjo'^ FOR 

Z=Zq and <I>(M)ocM-2-35 



Mup 

Mq 


VO 


^0 
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0.045 


0.091 


20 
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0.139 
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0.093 


0.209 
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0.121 


0.262 
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0.169 


0.354 


30 


0.186 


0.200 


0.424 
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0.246 


0.544 
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0.278 


0.635 


40 




0.359 


0.892 


45 


0.393 


0.445 


1.130 


50 




0.527 


1.347 


60 


0.598 


0.674 


1.739 



^ r/0 is taken directly from lVaccal 
1 119941) . we calculate 7}^ using 
the models presented in IVaccal 
jl994r , and we calculate r]^ 
based on the model presented in 
IMartins et al.l ||2005I1 . In our cal- 
culations, we cannot accurately 
report t/q for upper masses above 
60 Mq due to the lack of data 
above this mass in IMartins et all 

unni). 



Hibbard, J. E., van der Hulst, J. M., Barnes, J. E., & Rich, R. M. 

2001, AJ, 122, 2969 
Hills, J. G. 1980, ApJ, 235, 986 

Hoogerwerf, R., dc Bruijne, J. H. J., & de Zeeuw, P. T. 2001, A&A, 
365, 49 

Hunter, D. A., Elmegreen, B. G., Dupuy, T. J., & Mortonson, M. 

2003, AJ, 126, 1836 
Ivanov, G. R. 1996, A&A, 305, 708 
Kennicutt, Jr., R. C. 1989, ApJ, 344, 685 
— . 1998, ARA&A, 36, 189 

Kennicutt, Jr., R. C., Calzetti, D., Walter, F., Helou, G., 
Hollenbach, D. J., Armus, L., Bendo, G., Dale, D. A., Draine, 
B. T., Engelbracht, C. W., Gordon, K. D., Prescott, M. K. M., 
Regan, M. W., Tliornley, M. D., Bot, C., Brinks, E., de Blok, E., 
de Mello, D., Meyer, M., Moustakas, J., Murphy, E. J., Sheth, 
K., & Smith, J. D. T. 2007, ArXiv e-prints, 708 

Konstantopoulos, I. S., Bastien, N., Smith, L. J., & Westmoquette, 
M. S. 2007, in preparation 

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 

Lata, S., Pandey, A. K., Sagar, R., & Mohan, V. 2002, A&A, 388, 
158 

Leitherer, C., Schaerer, D., Goldader, J. D., Delgado, R. M. G., 
Robert, C., Kune, D. P., de Mello, D. P., Devost, D., & Heckman, 
T. M. 1999, ApJS, 123, 3 

Lejeune, T., & Schaerer, D. 2001, VizieR Online Data Catalog, 
6102, 

Leonard, P. J. T., & Duncan, M. J. 1990, AJ, 99, 608 

Martins, P., Schaerer, D., & HiUier, D. J. 2005, A&A, 436, 1049 

Mendes de Oliveira, C., Cypriano, E. S., Sodre, Jr., L., & 
Balkowski, C. 2004, ApJ, 605, L17 

Meurer, G. R., Hanish, D. J., Perguson, H. C., Knezek, P. M., 
Kilborn, V. A., Putman, M. E., Smith, R. C., Koribalski, B., 
Meyer, M., Oey, M. S., Ryan-Weber, E. V., Zwaan, M. A., 
Heckman, T. M., Kennicutt, Jr., R. C., Lee, J. C., Webster, 
R. L., Bland-Hawthorn, J., Dopita, M. A., Freeman, K. C., 
Doyle, M. T., Drinkwater, M. J., Staveley-Smith, L., & Werk, 
J. 2006, ApJS, 165, 307 

Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64 

Meynet, G., & Maeder, A. 2000, A&A, 361, 101 

O'Donnell, J. E. 1994, ApJ, 422, 158 

Oosterloo, T., Morganti, R., Sadler, E. M., Ferguson, A., van der 
Hulst, T., & Jerjen, H. 2004, in lAU Symposium, ed. P.-A. Due, 
J. Braine, & E. Brinks, 486—1- 

Oosterloo, T., & van Gorkom, J. 2005, A&A, 437, L19 



Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and 
active galactic nuclei (Research supported by the University 
of California, John Simon Guggenheim Memorial Foundation, 
University of Minnesota, et al. Mill Valley, CA, University 
Science Books, 1989, 422 p.) 

Rafelski, M., & Zaritsky, D. 2005, AJ, 129, 2701 

Ryan-Weber, E., Webster, R., & Bekki, K. 2003a, in Astrophysics 
and Space Science Library, Vol. 281, The IGM/Galaxy 
Connection. The Distribution of Baryons at z=0, ed. J. L. 
Rosenberg & M. E. Putman, 223-+ 

Ryan-Weber, E. V., Meurer, G. R., Freeman, K. C, Putman, M. E., 
Webster, R. L., Drinkwater, M. J., Perguson, H. C, Hanish, D., 
Heckman, T. M., Kennicutt, Jr., R. C, Kilborn, V. A., Knezek, 
P. M., Koribalski, B. S., Meyer, M. J., Oey, M. S., Smith, R. C., 
Staveiey-Smith, L., & Zwaan, M. A. 2004, AJ, 127, 1431 

Ryan-Weber, E. V., Webster, R. L., & Staveley-Smith, L. 2003b, 
MNRAS, 343, 1195 

Sakai, S., Kennicutt, Jr., R. C, van der Hulst, J. M., & Moss, C. 
2002, ApJ, 578, 842 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Schmidt-Kaler, T. 1982, in Landolt-Bornstein, New Series, Group 
VI, Vol. 2, ed. K. Shaifers & H. H. Voigt, 1— |- 

Sirianni, M., Jee, M. J., Bem'tez, N., Blakeslee, J. P., Martel, A. R., 
Meurer, G., Clampin, M., De Marchi, G., Ford, H. C, Gilliland, 
R., Hartig, G. F., lUingworth, G. D., Mack, J., & McCann, W. J. 
2005, PASP, 117, 1049 

Skillman, E. D. 1987, in Star Formation in Galaxies, 263-266 

Teplitz, H. I., Siana, B., Brown, T. M., Chary, R., Colbert, J. W., 
Conselice, C. J., de Mello, D. P., Dickinson, M., Ferguson, H. C, 
Gardner, J. P., & Menanteau, F. 2006, AJ, 132, 853 

Thilker, D. A., Bianchi, L., Meurer, G., Gil de Paz, A., Boissier, S., 
Madore, B. P., BoseUi, A., Ferguson, A. M. N., Mufioz-Mateos, 
J. C, Madsen, G. J., Hameed, S., Overzier, R. A., Forster, 
K., Friedman, P. G., Martin, D. C, Morrissey, P., Neff, S. G., 
Schiminovich, D., Seibert, M., Small, T., Wyder, T. K., Donas, 
J., Heckman, T. M., Lee, Y.-W., Milliard, B., Rich, R. M., Szalay, 
A. S., Welsh, B. Y., & Yi, S. K. 2007, ApJS, 173, 538 

Vacca, W. D. 1994, ApJ, 421, 140 

Vacca, W. D., Garmany, C. D., & ShuU, J. M. 1996, ApJ, 460, 914 
Walter, P., Martin, C. L., & Ott, J. 2006, AJ, 132, 2289 
Whitmore, B. C, Zhang, Q., Leitherer, C, Pall, S. M., Schweizer, 
F., & Miller, B. W. 1999, AJ, 118, 1551 



